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Abstract. We solve the dynamics of an open quantum system where N strongly driven two-level atoms are 
equally coupled on resonance to a dissipative cavity mode. Analytical results are derived on decoherence, 
entanglement, purity, atomic correlations and cavity field mean photon number. We predict decoherence- 
free subspaces for the whole system and the N-qubit subsystem, the monitoring of quantum coherence 
and purity decay by atomic populations measurements, the conditional generation of atomic multi-partite 
entangled states and of cavity cat-like states. We show that the dynamics of atoms prepared in states 
invariant under permutation of any two components remains restricted within the subspace spanned by 
the completely symmetric Dicke states. We discuss examples and applications in the cases N = 3,4. 

PACS. 42.50.Pq Cavity quantum electrodynamics; micromasers - 03.65.Yz Decoherence; open systems; 
quantum-statistical - 03.67.Bg Entanglement production and manipulation methods 



1 Introduction 

Cavity quantum electrodynamics (QED) concerns the in- 
teraction of atoms (or ions) with a quantized radiation 
field in a microwave or optical cavity [UIUGI]. The basic 
physical principles are quite well understood [4] and ac- 
curately tested [U[6,7J. Impressive advances in the experi- 
mental control on these systems allow to study fundamen- 
tal issues in quantum physics such as entanglement and 
decoherence, related to the nonlocal correlations of com- 
posite systems at the microscopic level and to the bound- 
ary between the quantum and the classical descriptions. 
These issues have raised a huge interest due to the poten- 
tialities of the peculiar quantum behavior for applications 
in quantum information processing, communication, and 
computation [|8j and the necessity to protect quantum co- 
herence from noisy environments [9lll0j . In this framework 
it is important to investigate nontrivial solvable models, 
representing somewhat idealized versions of systems im- 
plementable with the present cavity QED technology. In 
this paper we exactly solve the dynamics of an open multi- 
partite system, where N two-level atoms are equally cou- 
pled on resonance to a dissipative cavity mode and coher- 
ently driven by a strong external field. Cooled, trapped, 
deterministically loaded atoms pdl[T2iri3] and trapped 
ion systems [14 j ri5 ] in optical cavities, as well as Rydberg 
atoms crossing microwave cavities [T], appear as the most 
promising candidates for implementations. 
We present a quite compact solution of the open system 
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dynamics derived by phase-space techniques [IB], used in 
previous works [17, 18 . The N-atom subsystem can be de- 
scribed as a pseudo-spin system where the independent 
coupling of each atom to the cavity combines with the in- 
variance of dynamics under permutation of any two atoms. 
A peculiarity of the system is that this description does 
not hold in the standard (energy or computational) basis, 
but in a rotated one. The permutational invariance reflects 
in the atomic coupling to the environment, leading to 
the existence of both global and atomic decoherence-free 
subspaces (DFS) 19,20]. In the latter case an initial N- 
qubit entanglement remains protected and available e.g. 
for quantum memories or quantum processors [2TU22] . The 
structure of the general solution allows predicting a way 
to monitor the decay of quantum coherence and purity 
by measurements of atomic probabilities. In the limit of 
unitary dynamics these measurements can conditionally 
generate mesoscopic cat-like states of the cavity field. The 
preparation of atoms in states which are invariant under 
the exchange of any atom pair restricts the atomic dynam- 
ics to the subspace spanned by the completely symmetric 
Dicke states |23j , some of which are genuinely multipartite 
entangled states [53]. This further simplifies the descrip- 
tion of system dynamics in the important cases of atoms 
all prepared in a same state. Selected results for three and 
four atoms provide further insight on system and sub- 
systems dynamics, including preservation and conditional 
generation of multipartite entanglement. 
In Sect. [5] we introduce the model. The analytical solution 
of system dynamics is derived Sect. |3] where applications 
to N — 3, 4 are also reported. Results for transient and 
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steady state regimes are discussed in Sect. |4j and multi- 
partite entanglement protection is the object of Sect. El 
The main results are summarized in the Conclusions. 



2 The physical model 

We consider a set of N two-level atoms interacting with 
a dissipative cavity field mode. The transition frequency 
u> a between excited and ground states, \e)i and \g)i (I — 
1, N), is the same for all the N atoms. A coherent exter- 
nal field of frequency uj a simultaneously drives the atoms 
during the interaction with the cavity mode of frequency 
uif [2"5]l26]. This kind of system is feasible in cavity QED 
experiments with two-level Rydberg atoms in a microwave 
cavity [1] or with three-level atoms effectively reduced 
to two levels interacting with an optical cavity [17], due 
to relevant advances recently achieved in cooling, trap- 
ping and deterministically loading atoms in optical cavi- 
ties [TT] . In both regimes atomic decays can be neglected, 
as we shall assume from now on. Similar dynamics could 
be also implemented by trapped ions interacting with a 
cavity mode [14] . 
The whole system Hamiltonian is 



N 



1=1 

i nf n — iuj a t ~t i „iw a t ±. \ 

+ U{e <j\ + e o\) , 



(1) 



where Q is the Rabi frequency associated with the coher- 
ent driving field amplitude, g the atom-cavity mode cou- 
pling constant taken equal for all atoms, a (cr) the field 
annihilation (creation) operator, 07 = |<?)z(e| (a\ = \e)i{g\) 
the atomic lowering (raising) operator, and & z j — |e)/(e| — 
I <?)/((? I the inversion operator. 

In the perspective of experimental implementation of our 
scheme we add the effects of cavity mode dissipation, while 
we focus on resonance conditions (uj a = luj) in order to 
derive an analytical solution of the system dynamics. 
Therefore, we must solve the following master equation 
(ME) for the statistical density operator p' N of the whole 
system 



Pn 



(2) 



where 



where the Hamiltonian ([I]) has been replaced by the time- 
independent Hamiltonian H 1 = Ho + Hi with 



Ho 
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(5) 



Now we consider the unitary transformation U(t) = efi Hot 
and we derive for the density operator p^ = Up l N tv the 
following ME: 



Pn 
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(6) 



where the transformed Hamiltonian H[ — UHiW can be 
written as 



JV 



#i(*) = ^a£[a 



- 2UU )cn + (1 + e 2lQt )a\ 
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h.c. 

(7) 



In the strong-driving regime for the interaction between 
the atoms and the external coherent field, fl 3> g, we 
can use the rotating-wavc approximation obtaining the 
effective Hamiltonian [25II26] 
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(8) 



We notice the presence of Jaynes-Cummings (aja + &ja)) 

as well as anti- Jaynes-Cummings (crt^t _|_ ajh) coupling 
terms of each coherently driven atom with the cavity field. 
Hereinafter we shall solve the master equation for the 
whole system density operator f>N(t) 



Pn 



-[^off, Pn] + C-fPN- 



(9) 



3 Analytical solution and system dynamics 

In order to solve the general TV-atom problem we intro- 
duce the collective atomic operator S x — | X^=i &x,l = 

2 YliLii^i +o"{)> so * na ^ ^ ne effective Hamiltonan assumes 
the simple form 



£fPN = o ftap'r 
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(10) 



is the standard Liouville superoperator which describes 
the dissipative decay of the cavity field mode, with the 
rate k, due to the coupling to a thermal bath at zero tem- 
perature. 

In the interaction picture the dissipative terms remain un- 
changed and the ME ([2]) can be rewritten as 



p I N = -^[H I ,p I N } + C f p I N 



(4) 



We recall that the eigenstates of the spin operator a x ,i are 



the rotated states |± 



_ \gh±\eh 



i = ""-^ " where a Xt i\±)i = Af\^)i 
with Af = ±1. For the whole atomic subspace we consider 
the basis of 2^ states {|i)jv} where any \i)n is an eigen- 
statc of the collective spin operator S x . The correspond- 
ing eigenvalue Sj = (1/2) YI1L1 1S na ff tne difference 
between the number of |+) and |— ) components of state 
\i)ff, regardless of the exchange of any qubit pair, and it 
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can assume N + 1 values from —N/2 to N/2 with steps 
\Asi\ — 1. We notice that the eigenvalues have a de- 
generacy order given by n(si) = ( N / 2+Si )\(N/2- Si y. tnat is 
greater than one if -N/2 < s, < N/2. 
The general solution of the ME ([9]) can be derived by intro- 
ducing the decomposition of the density operator /3_/v(i) = 

2 N 

Tlii j=i Jv(i|/5jv(*)|j)jv|i)jv(j| on the N-atom rotated ba- 
sis |z)jv, so that it is equivalent to the following set of 
2 2N uncoupled evolution equations for the field operators 
pN,ij = N(i\pN(t)\j) N 



PN,ij = -ig[si(a + a r )p N4j - Sjp N ,ij(a + a')] 

+ £fpN,ij- 



(11) 



Equation can be solved by a combination of phase 
space techniques [16] with the method of characteristics 
[27 . Starting from the cavity in the vacuum state |0) and 
the atoms in any pure state 



1^(0)^ = 10}®^ 



CN,i\VN 



(12) 



the initial states of eq. (|T2|) containing only the corre- 
sponding n(0) = A!/[(A/2)!] 2 atomic eigenstates |i)jv- 
Let us consider for example the case of N = 4 atoms. 
The DFS is spanned by the tensor product of the cavity 
vacuum state and n(0) = 6 states with the same num- 
ber of |+) and |— ) components (see Tabic [T]). It preserves 
any initial global state within this subspace, protecting 
any entangled atomic preparation for quantum informa- 
tion purposes. 

Another interesting feature of the case with even N fol- 
lows from the presence in eq. (I13[) of terms with the cav- 
ity field in the vacuum state. Namely, if the optical cavity 
field is accessible to measurements, the absence of a re- 
sponse by an on/off detector generates a pure N-qubit 
state, that can be a multipartite entangled state. We con- 
sider again the N — 4 atoms case. Starting e.g. from 
the four atoms prepared in the ground state, the den- 
sity operator of eq. (I13p contains a timc-indcpcndcnt part 
(V6/4)(|0}<8|#}„), where \W) a = (1/V6)(| + + --) + | + 

-+-> + !+— +) + |-++-> + |-+-+> + |—++». 

Hence a null measurement of the optical cavity field gener- 
ates the pure 4-qubit state \\t r ) a whose entanglement prop- 
erties will be discussed later. 



with the normalization condition I Cjv,i | = 1, we ob- 

tain the compact solution for the whole system density 
operator 



3.1 Completely symmetric Dicke states 



M*) = E wU/(«)] (Si - S3 " )2 x (13) 

i,j=l ' 

x \-2s z a(t))(-2s 3 a(t)\ ® [*>jvO'I- 

We see that the dynamics correlates the eigenstates of S x 
with cavity field coherent states of amplitude proportional 
to 

o(t) = i|(l- e -i*). (14) 
The one-atom decoherence function 



/(t) = /x(t)e 2 l««l 2 



-t+'- 



.2|a(t)| 2 



(15) 



naturally splits into two parts: fi(t), which will appear in 

the atomic subsystem dynamics, and e 2 '"^*-" , that is the 
field states normalization. It is responsible for the decay of 
coherences and depends on the dimensionless parameters 
(g/k) 2 and kt. 

In order to evaluate the degree of mixedness of the state 
pjsi(t) we derive from eq. (fl~3| the purity 



m (t)=Tr[p 2 N (t)} 



2" 

E 



CN,i\ 2 \cN tj \ 2 [f(t)] 



2(s,-s j f 



(16) 



A remarkable consequence of the general solution of eq. JT3 
is the existence of a global DFS for any even value of the 
number N of atoms, when the eigenvalue Sj can assume 
the value zero. In this case there is no time evolution for 



In order to exploit the system dynamical invariance un- 
der exchange of any atom pair, we consider initial states 
(eq. (|12p ) having in the atomic part only symmetric states 
or symmetrized combinations of states |i)jy. In this case 
the atomic part of /3jv(i) remains confined in the subspace 
spanned by only N +1 (instead of 2^) states that we de- 
note as |y,s), where — ^ < s < ^ with steps \As\ = 1. 
We notice that the above states are analogous to the so- 
called completely symmetric Dicke states (CSD) in [23], 
written in the energy basis instead of the rotated one. For 
instance for N = 4 the states \2,s) with —2 < s < 2 
are the symmetrized combinations of the states listed in 
Table 1. All the previous treatment can be adapted corre- 
spondingly. In particular, starting from any superposition 
of CSD states 



N/2 N 

M0))n = |o) ® Yl ^ly> s > 

s=-N/2 



N/2 



with the normalization condition V* 
the general solution (fT3|) can be rewritten as 



-N/2 



J N,s\ 



N/2 

p N (t)= Yl b N,sb N y[f(t)} 
s,s'=-N/2 



(s-s'f 



(17) 



= 1, 



(18) 



N N 

x | - 2sa(t))(-2s'a(t)\ ® \-,s)(-,s'\. 

In this case the interaction correlates cavity field coherent 
states with atomic CSD states. These results include the 
important case of all atoms prepared in the ground state, 
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Table 1. Atomic DFSs for TV = 4 qubits. 

s ; n(si) |i) 4 
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Fig. 1. Time evolution of the purity of the whole system for 
different values of N and for the fixed dimensionless parameter 
g/k = 5. 



where 6jv )S = (l/v2^)u(s). In another relevant case, with 
all atoms prepared in the excited state, the only change 
in (fl"5|) is the replacement f{t) — > —f(t). In such cases the 
purity (|16p reduces to 



N/2 



Pn( 



(*) = ^f E /(*) 



(19) 



3.2 Subsystem dynamics 

We derive some general results on the cavity mode and 
the atomic subsystems, providing the time-dependent ex- 
pressions of the corresponding reduced density operators 
and purities, as well as of the mean number of photons in 
the cavity. 

If we trace the whole system density operator of eq. (|T3| 
over the atomic variables, we find the expression for the 
reduced density operator of the cavity field 



Pn 



f (t) = M 2 \ - 2sMt))(~2sMt)\ (21) 



i=l 



that is a statistical mixture of coherent states, and whose 
purity is 



VNj(t) = Tr[p 2 Nf {t)\ = 



E \c N4 f\c Ntj \ 2 [e^Y Si - s ^ 



(22) 



s,s'=-N/2 



From the density operator pN,f(t) we can derive the ex- 
pression of the mean photon number 



whose asymptotic value can be written in a closed form in 
terms of the gamma function r 



ita>(t) = Tr f \p N , f (t)tfa] = 4|a(t)| 2 £ a 



2|„ |2 
CN,i\ ■ 



(23) 



N/2 



N 



^ ~ 22N E ~ 2?N 

=-iV/2 



=0 



r(N- 



(20) 



In the case of all atoms in the ground state, cna = 1/ v^^, 
one obtains 



where I = s + N/2 and ( , ) is the binomial coefficient. In 
Fig. [1] we show the time evolution of the system purity 
(|19[) for a fixed value of g/k and different qubit numbers 
N = 1, 4, where the steady state values are pf 5 = 1/2, 
pSS = 3 / 8j ^SS = 5 / 16) M |S = 35/X28. The greater the 

value of N, the faster the decay of the global coherences. 
Varying the ratio g/k instead of N, the asymptotic behav- 
ior does not change whereas the decay is faster (slower) 
for increasing (decreasing) values of g/k. 
We remark that the atomic preparation in one of the CSD 
states is equivalent to the qubit encoding in the corre- 
sponding DFS. In particular, for even values of N, the 
CSD atomic state with s — and the cavity in the vac- 
uum state belong to a global DFS. 



a>(t) = 



>N-2 



1=0 



= «|o(t)| 2 (24) 



showing that each atom gives the same average contribu- 
tion to the cavity field. 

By tracing the whole system density operator over the field 
variables, we obtain the reduced atomic density operator 

PN,a(t) = E C N,iC* Ntj [A(t)}^- S ^\i) N (j\. (25) 

We notice that if the atoms are prepared in any super- 
position of eigenstates corresponding to a degenerate 
eigenvalue —N/2 < Si < N/2, the state does not evolve. 
Therefore we identify N — 1 atomic DFSs with dimension 
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Fig. 2. Comparison between the iV-qubit purity /ijv, a (i) in 
<I26|) evaluated at steady state and 1/2 that is the purity of 
a maximally mixed state. 



that is an atomic correlation (bunching) effect, whereas 
the joint probability of the other two outcomes is 3/16, 
showing an antibunching effect. By exploiting the above 
expressions (|2"T1) it is possible to monitor the two-atom 
decoherence function ff(t) measuring the sums P eee {t) + 
Pggg (t) or P eeg (t) + P egg (t) . Remarkably the N-qubit de- 
coherence originates from the one-atom decoherence func- 
tion [17j which can be monitored by atomic population 
measurements via the relation fi(t) — P g (t) — P e {t), as 
well as the N-qubit purity according to (|26jl . In that case 
the atom and the field can approach maximally entangled 
states in the limits kt 1 and (g/k) 2 ^ 1, and the en- 
tanglement (measured by the Von Neumann subsystem 
entropy) is also described by fi{t). 



n(si) greater than one, where an initial entanglement can 
be protected. The purity of state (|25|) is 

A*iV,a(t) = Tr[p 2 Na {t)\ = 

= E i^i 2 i^-i 2 [/i(t)] 2( -"^ 2 . 12,11 

The decay of atomic purity ruled by fi(t) is faster than 
the global purity decay (fTB")) . ruled by f(t). However the 
asymptotic behavior is the same and we can use the re- 
sult (f20|) in order to make a comparison with the case of 
maximally mixed states, whose purity is equal to 1/2 N . 
In Fig. [2] we see that the state is maximally mixed only 
for N = 1, where actually the atom becomes maximally 
entangled with cavity field [T7] . For any N > 1 the state is 
never maximally mixed due to the survival of coherences 
in the DFSs. Also we notice that the field purity (|22l) 
remains slightly larger than the atomic one because the 
decoherence function fi (t) is replaced by a non- vanishing 
exponential function. 

As an application of the atomic subsystem dynamics we 
rewrite the atomic density matrix eq. (|25p in the stan- 
dard basis for the case N — 3. Starting, for instance, from 
the three atoms in the ground state the diagonal matrix 
elements provide the following joint probabilities for the 
atomic level populations 



Peee (^) 
Peeg (^) 
Pegg (t) 

Pggg W 



1 

32 
1 

32 
1 

32 
1 

32 



10 -15/!(t)+6/?(i) -/?(*) 
2-f l {t)-2ft{t)+f 1 {t) 

2+h{t)-2ft{t)-m 

10 + 15/ 1 (i) + 6/ 1 4 (<) + fl{t) 



(27a) 
(27b) 
(27c) 
(27d) 



where eqs. (|27b[) and (|27c[) represent one third of the 
probability to detect, respectively, two atoms in the ex- 
cited state or in the ground state, independently from 
the atomic ordering. The three-atom probabilities (|2T[) are 
shown in Fig. [3] 

We can see that at steady state the joint probability 
that three atoms are in the same state is equal to 5/16, 



4 Transient and steady state results 

In the Hamiltonian limit, kt <ti 1, of small cavity decay 
rate and/or short interaction times, f(t) — » 1, a(t) — ► 
a(t) = if and p N (t) -> \§(t)) N {W(t)\ where the global 
cat-like state 



2 



\V(t))N = > v c Nti \ - 2si&(t)) ® \i) N . 



i=l 



(28) 



As an example, for N — 3 we consider the generation of 
the pure state |<^(i))3. Starting from three atoms in the 
ground state, so that c^i = 1/V&, we obtain for H< 1 
an evolved state that we rewrite in the standard atomic 
basis 



a) + 3\a) - |3a)) <g> |eee) 



1 

71 


(l"3a) 


3a) 


F3|-a 


3a) 


-\-a) 


3a) 


+ \-&) 



- 3|a) + |3a» ® \ggg) 
\a) + 1 35)) ® {\eeg) + 



|ege) + \gee)) 



(29) 

where for brevity we have defined a = a{t). We notice a 
superposition of mesoscopic cat-like states of the cavity 
field correlated with atomic states with the same number 
of ground (or excited) atoms, which are two fully separable 
and two entangled 3-qubit states (a W and an inverted- 
W state) [2S]. An interesting consequence of eq. (f2"9"| is 
that a simultaneous detection of the three atoms in any 
state prepares the cavity field in the corresponding cat- 
like state. In Fig. 2] we show the Wigncr function that 
describes in phase space the cat-like state generated for 
atomic detections in the ground state. 
After the transient the coupling of the field to the envi- 
ronment introduces in the solution (|13j) the field-atoms 
coherences f(t),f 4 (t),...,f N (t). Note that these powers 
of the decoherence function can be obtained by the sub- 
stitution g —> Ng, which exactly reflects the independent 
interaction of each atom with the cavity field. 
In the steady state limit kt 3> 1 the density operator /3jv(i) 
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becomes a statistical mixture of the pure states superim- 
posed in the global cat-like state (f2"5)) generated in the 
transient 



Pn 



\cnA 2 \ - 2 Sl a ss ){-2 Sl a ss \ ® \i) N {i\, (30) 



where a ss = «f ■ The system (subsystem) purity at steady 
state was discussed in the previous section. 




Fig. 4. Wigner function W ggg of the cavity field state condi- 
tioned to the detection of the three atoms in the ground state, 
for parameters values kt — 0.05 and g/k — 110. 



5 Protection of multipartite entanglement 

Now we recall some concepts and tools in order to an- 
alyze the multipartite entanglement properties of some 
atomic states encoded in DFSs. The 3-tangle measure in- 
troduced in [29 evaluates the amount of entanglement 



shared by all the three qubits through the quantity T123 = 
C 2 2 + C13 — ^1(23)' wnere C-ij is the concurrence of the 
qubit pair (i, j). A generalization to the case of N qubits 
(with N even) was given in [3D] by the iV-tangle measure 
defined as t n = |(V#)| 2 with |$) = af N \ip*), where \tp) is 
the generic iV-qubit state, its complex conjugate and 
(j y one of the Pauli matrices. Another useful tool is the 
residual bipartite entanglement measure (see [55j) which 
evaluates the robustness of entanglement against the loss 
of information; this measure is provided, for instance, by 
the average squared concurrence C 2 calculated for any two 
residual qubits when the other N — 2 are traced out. 
Following the previous analysis concerning the DFSs we 
first consider the atomic subsystem for TV = 3. In this 
case we find that the three qubits do not evolve in time if 
they are prepared in any of the four decoherence-free CSD 
states 



|3/2,3/2) = 
|3/2,l/2) = 

|3/2,-l/2) = 
|3/2,-3/2)r 



l + +-) + l + -+> + |-++> 



V3 



+- + 



V3 



(31a) 
(31b) 

(31c) 
(31d) 



Two of them, (|31ap and ( 31d|) . are manifestly separable. 
The other two states, (I31b[) and (|31cp . show interesting 
entanglement properties. They have no full tripartite en- 
tanglement (T123 = 0) according to the 3-tangle mea- 
sure. However each qubit pair retains the maximal resid- 
ual bipartite entanglement C 2 = 4/9. These kind of states 
show a multipartite entanglement characteristic of W-like 
states. 

Let us now investigate the dynamics of four qubits which 
presents five decoherence-free CSD states, including two 
separable states |2,±2), and the multipartite entangled 
states |2,±1) and |2,0). The relevance for applications in 
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quantum information processing is that the state |2, 0) 
turns out to be maximally entangled according to the 
4-tangle measure (r 4 = 1), whereas the states |2, ±1) 
have no four-partite entanglement (t 4 — 0), but each of 
them exhibits an equal maximal reduced bipartite entan- 
glement, C 2 — 1/4 (W-like states), by tracing over any 
qubit pair. We remark that all CSD states of the type 
|iV/2, s) with s = ±(N — 2)/2 have entanglement proper- 
ties similar to that of states \Wn) introduced in [28] . By 
tracing the atomic density operators I-/V/2, s)(N/2, s\ over 
any N — 2 parties we always obtain the reduced density 
operators for the bipartite system 

p ± = A (2|#-><*-| + (N - 2)| ± ±)<± ± |) . (32) 

Hence for the average squared concurrence we simply ob- 
tain the value CP = (2/N) 2 . 

We further notice that we can rewrite the states |2,±1) 
and |2,±0) as 

|2,±1) = ± ±)i 2 |#->34 + |* _ )ia| ± ±> 34 ) 

|2,0) = -L(|<?-) 12 |<?-) 34 + |#-) 13 |#-)a4+ (33) 

+ |*~)i4|#~)a3) 

thus generalizing the results derived in [18j where we showed 
that the two-atom maximally entangled Bell state = 
775 (H — )ij + I — 0*j) an d the two separable states | ± ±) 
do not evolve in time. 

Another interesting application is to encode the four qubits 
in some states of a special basis called Bell gem [51] . 
which is a generalization of the well known Bell basis 
= (1/V2)(|0 5 >±|ec» and |lP±) = (1/V2)(\ge)±\eg)). 
It is composed by maximally entangled states, according 
to the 4-tangle measure (t 4 = 1), which can be obtained by 
simple quantum logic circuits starting from four unentan- 
gled qubits in the computational basis. Let us consider the 
cavity field prepared in the vacuum state |0) and the four 
qubits in one of the last three elements of the Bell gem 
(l/V2)(|<P+#+) - (l/\/2)(|#-$-) ± I* - !? - )). 

The whole system does not evolve in time because these 
three initial atomic states belong to the DFS correspond- 
ing to the eigenvalue Sj = 0, thus maintaining the maxi- 
mum multipartite entanglement in the atomic subsystem. 

6 Conclusions 

We have solved the dynamics of a feasible cavity QED sys- 
tem where N strongly driven two-level atoms are equally 
and resonantly coupled to an optical field mode in contact 
with an environment and initially in the vacuum state. 
For negligible atomic decay we have derived a compact 
solution of the open system master equation in terms 
of coherent field states, atomic pseudo-spin states, and 
suitable decoherence functions. We have derived and dis- 
cussed a number of exact results on system and subsys- 
tems dynamics which are also of interest for quantum 



information applications, including decoherence-free sub- 
spaces and multipartite (iV-qubit) entanglement protec- 
tion. In addition we have suggested a way to monitor de- 
coherence by atomic population measurements. For atoms 
prepared in symmetric states with respect to the exchange 
of any atom pair, including the physically important prepa- 
ration in the same (ground or excited) state, the dynamics 
is entirely expressed in terms of symmetric Dicke states. 
Applications in the cases with N = 3,4 have been dis- 
cussed. 
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